Quantum tunneling of the Neel vector in antiferromagnetic [3 x 3] grid molecules 

0. WaldmanrE 

Department of Chemistry and Biochemistry, University of Bern, CH-3012 Bern, Switzerland 

(Dated: February 2, 2008) 

Based on numerical calculations it is shown that the antiferromagnetic grid molecule Mn-[3 x 3] 
is a very promising candidate to experimentally detect the phenomenon of quantum tunneling of 
the Neel vector. 
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Molecular nanomagnets such as Mni2 and Fes have 
attracted much interest as they show quantum tunneling 
of the magnetization at low temperaturesi For antifer- 
romagnetic (AF) molecular wheels, a class of magnetic 
molecules in which an even number of ions is arranged 
in a ringlike fashion, 2 a different scenario associated with 
quantum tunneling of the Neel vector (QTNV) was pre- 
dicted theoretically, 3 but so far could not be observed ex- 
perimentally. The identification of a magnetic molecule 
which enables the detection of QTNV would represent 
a major achievement, of interest also for applications in 
quantum technologies. 

The experimental task to detect QTNV in molecular 
wheels is challenging: high sensitivity, low temperatures 
and strong magnetic fields are mandatory. A technique 
such as electron spin resonance (ESR) could meet these 
requirements, but does not couple to the Neel vector be- 
cause of symmetry reasons4£ One way out would be to 
have a " sensor spin" which is coupled to only one of the 
AF sublattices^ Its dynamics then will reflect the dy- 
namics of the Neel vector. In principle, the nuclear spin 
of a metal center or a nearby proton might act as a sen- 
sor spin. Another solution would be to replace one of 
the magnetic ions on the ring by a dopant ion with dif- 
ferent spin*! The resulting excess spin also reflects the 
dynamics of the Neel vector - with the advantage of be- 
ing detectable by ESR. Efforts in this direction are under 
way, but the synthesis of appropriate modified wheels is 
quite difficult. 

Recently, a supramolecular Mn-[3 x 3] grid was 
investigated^ Here, nine spin-5/2 Mn(II) ions occupy 
the positions of a regular 3x3 matrix, held in place by 
a lattice of organic ligands [see inset of Fig. [Ja)]. The 
unusual magnetic properties suggested that the [3 x 3] 
grid can be regarded as an AF "ring" doped with a cen- 
tral Mn ioni Indeed, the spin Hamiltonian of a [3 x 3] 
grid, which for an idealized structure reads^ 

H = -j r \Y, s > ' S *+i + s « ' Si J +d r J2 sL 

\i=l J i=l 

~Jc (Sa + S4 + Se + Ss) ■ Sg 
+D c Sl tg +gn B S-B, (1) 

consists of a part involving only the ring of the eight pe- 
ripheral metal ions, a part related to the central ion, and 
terms representing an interaction between these two sets 
of spin centers [Si is the spin operator of the ith ion with 



spin s = 5/2, spins at "corners" are numbered 1, 3, 5, 7, 
those at "edges" 2, 4, 6, 8, and the central spin is num- 
bered 9, see Fig.^a); z denotes the axis perpendicular to 
the grid plane] . For Mn- [3 x 3] , an AF intraring coupling 
( Jr < 0) and uniaxial easy-axis anisotropy (Dr < 0) was 
found. 7,8 Furthermore, the central spin is coupled to only 
one of the sublattices of the peripheral ring, namely to 
the edge spins. Thus, magnetically the Mn-[3 x 3] grid 
should correspond to an AF ring of eight spin-5/2 centers 
with a sensor spin naturally built in. This work demon- 
strates that Mn-[3 x 3] is indeed a promising molecule 
to detect QTNV experimentally. 

The dimension of the Hilbert space is huge for Mn- 
[3 x 3]. However, the low-energy sector of an AF 
ring (with a small, even number of sites N) is very 
well described by approximating the wave functions by 
\PaPbSaS b SM) with S A = S B = Vs/2£±2iil Here, S A 
(Sb) denotes the total spin of sublattice A (B), and Pa, 
Pb abbreviate intermediate quantum numbers (omitted 
in the following). Physically, this approach works well 
because the internal spin structure due to the dominant 
Heisenberg interaction is essentially classical^ Applying 
this approximation to the ring terms in Eq. Q , an effec- 
tive Hamiltonian for the Mn-[3 x 3] grid is obtained: 

H e ff = —JrSa ■ Ss + Dr{S 2 A z + S 2 Bz ) 

-J C S A ■ S 9 + D c Sl z + g^ B S ■ B. (2) 

Here, Jr = 0.526 J and Dr = 0A97D R . A (B) denotes 
the sublattice of corner (edge) spins. Figure ^ compares 
results calculated for H^ff and H. The agreement is 

quite good. H^ff is designed to reproduce the lowest- 
lying states optimally, higher-lying states are thus in- 
creasingly less well described. 

Remarkably, H^Jf works well in a wide range of 
Jc/\Jr\ values. In the strong coupling limit, \Jc\ 3> 
\Jr\, the wave functions should be approximated by 
\S b SqSbqSaSM), i.e., by coupling first the spins Sb 
and Sg, and then the resulting spin Sb9 with Sa- In 
the weak coupling limit, however, the appropriate wave 
functions would be \SaSbSabSqSAI), where first Sa and 
Sb are coupled together. As both wave functions lead to 
the same effective Hamiltonian, Eq. J2J), a broad range 
of Jc/\Jr\ is covered. This suggests a general strat- 
egy to construct an effective Hamiltonian: Sublattices 
as expected from the classical spin configuration are in- 
troduced, and then spins of each sublattice are replaced 
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FIG. 1: (a) Energies of the low- lying states with respect to the 
5 = 5/2 ground state as function of Jc for H^? (lines) and 
H (symbols) with Dr = Dc — in zero field. For H, only the 
energy gaps for the lowest states with S = 7/2, 9/2, 11/2 are 
displayed. The inset shows the structure of the Mn-[3 x 3] 
grid schematically, (b) Energy spectrum for H^jf (lines) and 
H (circles) as function of D = Dr = Dc for Jc ~ Jr at zero 
field. The spin multiplets get split by the anisotropy. The S = 
5/2 level for instance split in the sequence \M\ = 1/2, 3/2, 5/2. 
Twofold degenerate spin levels such as the S = 3/2 levels at 
about 3.5|Jfl| are not reproduced by H^f- They belong to 
an E band similar as in wheelsrS^ but are not relevant here 
as explained in the text. 



by the mean-field spin S S l = 1/Nsl J2i&SL S * ( n sl de- 
notes the number of spins of sublattice SL). The effect 
of weak quantum fluctuations may be accounted for by 
renormalizing the parameters of the effective Hamilto- 
nian, hence Jr and Dr in Eq. JSJ). 5 

Interestingly, a similarity of the single-molecule mag- 
net Mn ^ 1 ! 13 and a [3 x 3] grid turns up: Mni2 may be 
regarded as an octanuclear ring of Mn(III) ions doped 
by a central Mn(IV) tetramcr. The above considerations 
then lead directly to an effective Hamiltonian, which, in 
view of the success of this approach for the molecular 
wheels and the [3 x 3] grid, is expected to describe also 
the relevant states of Mni2 well. 

With H^ f 3 , QTNV in Mil- [3 x 3] can now be readily 
analyzed numerically. The energy spectrum as function 
of magnetic field is shown in Fig. [3] for several field ori- 
entations and Jr = Jc — —5 K, Dr = Dc = —0.14 K 
(9 denotes the angle between magnetic field and z axis; 
the azimuthal angle tp has no effect, as is evident from H, 
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FIG. 2: Energy spectrum vs magnetic field for different orien- 
tations 9 of the magnetic field (the energy of the lowest state 
was set to zero at each field). For zero field, the states with 
energy below 3 K belong to the S — 5/2 ground state multi- 
plet (which exhibits a hard-axis zero- field splitting), those at 
« 12.5 K to the next higher S = 7/2 multiplet, and those at 
w 25 K to the S = 9/2 multiplet. 



H e f~f). These parameters are supported by recent inelas- 
tic neutron scattering experiments^ and torque mea- 
surements at very low temperatures. 8 The semiclassical 
tunneling action S /H = NS t y / 2D R J J (QTNV sets in for 
So/h 3> 1, see Ref. 0) is then estimated to So/h = 4.7, 
which is even larger than for Feio (So/h = 3.3) and CsFeg 
(S /h = 3.9)M This demonstrates that Mn-[3 x 3] is in- 
deed a favorable candidate with respect to QTNV. Some 
following points should be noted: 

1) Concerning the dependence on Jc/\Jr\, the states 
should be divided into the energetically lowest states for 
S = 5/2, 7/2, 9/2, . . ., and the remainder of the states (in 
the terminology of Refs.llOl and ll2l these states form the 
L band, while the remaining states belong to either the 
E band or the quasicontinuum) . Obviously, only the first 
set of states (the L band) is relevant for the low- lying part 
of the energy spectrum in magnetic fields. As soon as Jc 
assumes significant values, the relative energy spacings of 
these states, and therewith the relative field positions of 
the ground-state level crossings, are close to those found 
in the strong-coupling limit \Jc\ S> \Jr\, even if \Jc\ is 
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much smaller than \Jr\. The weak-coupling limit is real- 
ized only for \ Jc/Jr\ < 0.002. Thus, the results shown 
in this work for Jq = Jr are actually characteristic for 
the strong coupling limit. 

2) The energy gap between ground and first excited 
state, A(B), shows a | sm(ng^BB/ Jr + a)|-like oscilla- 
tory field dependence for high fields in the grid plane 
[Fig. 12(d)]. Such a behavior was identified as a charac- 
teristics of QTNV in the (modified) molecular wheels. 3 6 
a accounts for a shift in field. 

3) Finally, A(B) periodically drops to zero indicating 
level crossings (LCs) at fields B m (m <E N)£±i Such LCs 
trivially exist for 9 = 0°. However, for canted magnetic 
fields gaps appear at each LC field, i.e., the ground state 
LCs become anticrossings signaling level mixing (Fig.[2J). 
Starting with 9 = 0°, the gaps at the LC fields first in- 
crease, reach maximal values at around 9 = 45°, and then 
shrink to disappear again at 9 = 90°. In the (unmodi- 
fied) wheels, A(B m ) = is enforced by symmetry for all 
field directions This symmetry is absent in [3 x 3] 
grids, and for fields with 9^0° gaps open at the LCs 
due to the action of the uniaxial magnetic anisotropy: In 
a coordinate frame with z' parallel to the field, one ob- 
tains DS 2 = D'S' 2 + E'{S 12 - S' y 2 ) - G'{S' X S' Z + S' Z S' X ). 
The G' term mixes the two levels involved in a LC (for 
which | AM' | = 1), leading to the above behavior as 
G' = sin(0) cos(9)D /3. This angular behavior of A(B) 
has important implications for the QTNV scenario (vide 
infra). 

In order to analyze QTNV, the matrix elements be- 
tween the ground state |0) and the 20 next higher lying 
states |n) were calculated for the z component of the Neel 
vector |(0|iV z |n)|, the total spin |(0|iS^|n)|, and the cen- 
tral spin KOjSg^ln)! (n = 1,...,20 numbers the states 
by their energy at each field). The Neel vector was de- 
fined as N = S^4 — Sr + Sg. The results as function of 
magnetic field are shown in Figs. |3(a)-(c) for 9 = 90°. 

For the Neel vector, the excitation between ground 
and first excited state, |(0|iV z |l)|, which corresponds to 
QTNV, is dominant [Fig. ^a)]. The energy gap ALB) 
is thus identified as the tunneling gap due to QTNV 
(for fields beyond the first LC). Higher-lying excita- 
tions, in particular, to the second excited state, are also 
present. Accordingly, the correlation function vs time, 
{0\N z (t)N z (0)\0) oc J2 n \{0\N z \n)\ 2 exp(-jrE n t), shows a 
dominant oscillation at the tunneling frequency A/h dis- 
turbed by higher frequent oscillations. As the contribu- 
tion of higher-lying excitations is minimal for magnetic 
fields inbetween LCs [where the tunneling gap A(B) as- 
sumes maxima, see Fig. [2(d)], the dynamics of the Neel 
vector is best described by QTNV at these fields. 

In S z , three excitations have significant intensities 
[Fig. [3(b)] . The excitation with maximum at ss 1 T cor- 
responds to a transition within the S = 5/2 ground-state 
multiplet. The excitation gaining increasing intensity 
from a 5 T onwards corresponds to a transition to the 
Zeeman-split level, which is visible in the energy spec- 
trum at energies r* gfiBB [Fig. [2(d)]. And finally, the 
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FIG. 3: Absolute matrix elements between the ground and 
the 20 next higher- lying states for the z component of (a) 
the Neel vector, (b) the total spin, and (c) the central spin 
for fields in the plane of the grid. Arrows and dashed lines 
in panel (a) indicate the positions of LCs. (d) Matrix ele- 
ments for the Neel vector with slightly canted field. The inset 
magnifies the field range 11 - 13 T for the first three matrix 
elements. 



excitation with an approximate B~ 2 field dependence 
corresponds to the transition between the ground and 
first excited state reflecting exactly QTNV. 6 

The matrix elements for the central spin basically just 
behave as those of the Neel vector [Fig. [3(c)]. This is of 
course expected, but demonstrates unambiguously that 
the central spin of a [3 x 3] grid indeed acts as a sensor 
spin for the dynamics of the spins on the peripheral ring. 

The results resemble those for molecular wheels: The 
field dependencies of A(B), |(0|5 Z |1)|, |(0|V 2 |1)| (and 
|(0|Vz|2)|) show similar behavior as for the modified and 
unmodified wheels&ii As shown above, the strong cou- 



4 



pling limit is effectively realized in Mn-[3 x 3]. Then 
the sublattice spin Sg and the central spin Sg act as 
a combined larger spin Sb9, and H^j? simplifies to 
H' = -J>S A -S B9 + D' A S 2 A z + D' B S 2 mz +g f i B S-B, which 
is exactly the effective Hamiltonian of a modified wheel, 6 
Thus, magnetically Mn-[3 x 3] behaves like a modified 
AF wheel. 

The new feature observed here is that, along with the 
opening of gaps at the LCs, the behavior of the matrix 
elements changes drastically for 9 ^ 0,90°, Fig. Efd) 
(preliminary calculations confirmed similar behavior for 
the modified wheels). At the LC fields, |(0|7V Z |1)| (and 
|(0|iV 2 |2)|) vanishes and |(0|JV Z |3)| becomes the dominant 
matrix element. Apparently, the opening of gaps at the 
LCs and the rise of |(0|7V Z |3)| results in a breakdown of 
the QTNV scenario at these fields. Thus, the picture 
of QTNV needs an extension near LCs in the case of 
Mn-[3 x 3] or modified wheels. However, importantly, 
in those regions of the magnetic field where the tunnel- 
ing gap A(B) is large, QTNV persists even for sizeable 
canting angles. 

Finally, the prospects of an experimental observation 
of QTNV in Mn-[3 x 3] is discussed. Having shown that 
the central spin acts as sensor spin and the analogy with a 
modified wheel, the considerations of Refs.Uy are appli- 
cable to Mn-[3 x 3]. Only main points will be addressed 
here, further details may be found in these works. 

QTNV can be detected by the highly sensitive ESR 
technique, as is evident from Fig. Efb) noting that the 
ESR intensity is proportional to |(0|S 2 |1)| 2 for 9 = 90°. 
For Mn-[3 x 3], the most favorable field regime would be 
10-20 T, see Fig.EJa). In this regime, the ESR intensity 
is smaller by two orders of magnitude compared to the 
maximal intensity at B = 0, but is still large enough to 
be detectable with today's ESR spectrometers, especially 
as rather large high-quality single crystals are available. 15 
The magnetic field should be oriented close to the plane 
of the grid, but a canting of about 5° is well tolerable 



if QTNV is measured near the tops of the tunneling gap 
A(.B), i.e., at fields inbetween LCs. In a continuous-wave 
ESR experiment, the linewidth of the signal provides an 
upper limit for the decoherence rate T s, which allows one 
to test the coherence condition Ts < A/Ti., or if QTNV 
in Mn-[3 x 3] is coherent, respectively^. 

Nuclear magnetic resonance (NMR) experiments are 
an alternative. As shown in Ref. |J, QTNV leads to an 
absorption peak in the NMR signal at u — A/ft, but 
with an intensity significantly reduced as (A/ A) 2 (A is 
the hyperfine coupling constant). For molecular wheels 
it was suggested to measure near LCs, where A is small, 
in order to gain intensity. But for Mn- [3 x 3] , this would 
require an unrealistically accurate orientation of the mag- 
netic field perpendicular to the grid plane. Fortunately, 
the sensitivity of Mn 55 -NMR is 3xl0 5 times larger than 
of Fe 57 -NMRji& resulting in detectable signals also near 
the maxima of A(B) in Mn-[3 x 3]. And, as mentioned 
already, relatively large single crystals are available^ It 
is added that the NMR signal due to QTNV can be easily 
identified by its field dependence which differs markedly 
from that of the Lamor frequency. 

In conclusion, the molecular grid Mn-[3 x 3] is an ex- 
citing prospect for observing quantum tunneling of the 
Neel vector experimentally: It has a sensor spin natu- 
rally built in which enables a detection of tunneling by 
ESR, and it features the highest value of the tunneling 
action Sq/H known to date. The issue of coherence is 
of great importance for any potential applications. In 
view of tunneling gaps of 1 K, coherent QTNV might be 
realized in Mn-[3 x 3]. 
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